Flux Flow Instabilities in Microstructured Amorphous NboyGeos Thin Films 
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We report measurements of the electric field vs. current density [-E(J)] characteristics in the 
mixed state of amorphous Nbo.yGeo.a microbridges. Close to the transition temperature Tc the 
Larkin-Ovchinnikov theory of nonlinear flux flow and the related instability describes the data 
quantitatively up to ~ 70 % of the upper critical magnetic field Bc2 and over a wide electric field 
range. At lower temperatures the nonlinearities of E{J) can be described by electron heating which 
reduces Bc2 and leads to a second type of flux flow instability, as shown by a scaling analysis of the 
high-dissipation data. 
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PACS numbers: 74.78.Db, 74.40.+k, 74.25. Qt 

It was predicted by Larkin and Ovchinnikov (LO) that 
the E{J) curves of a "dirty" superconductor in the mixed 
state may exhibit a steep increase long before the depair- 
ing current density is reached»i This behaviour, called 
LO flux flow instability (FFI), appears when the energy 
of quasiparticles in driven vortex cores becomes large 
enough for their escape into the surrounding superfluid. 
As a consequence, the vortices shrink and the vortex mo- 
tion viscosity is reduced, resulting in an increase of the 
dissipation at the fixed J. According to LO, the above 
mechanism is strong at temperatures T close to Tc and 
the corresponding FFI occurs when E reaches a critical 
value Ei proportional to magnetic field B, which implies 
a i?- independent critical vortex velocity Ui = Ei/ B. An 
extension of the LO theory by Bezuglyj and Shklovskij 
(BS)^ who took into account heating effects due to a fi- 
nite rate of heat removal to the bath, limited the pure, 
nonthermal LO FFI to a i? considerably smaller than 
Bc2- In order to explain a i?-dependent Ui, other modifi- 
cations of the LO theory explored nonlocality of the spa- 
tial distribution of the excitations, caused by the quasi- 
particle energy relaxation length being smaller than the 
intervortex spacingnSi Explanations of the FFI beyond the 
original or modified LO picture were sought in dynamic 
vortex lattice crystallisation,'* depinning phenomenajS, 
appearance of hot spots,*' and recently in vortex core ex- 
pansion due to electron heating at low temperatures ii 
Irrespective of the microscopic origin of the FFI, its dis- 
tinct feature is an E[J) region just above Ei where theory 
predicts dE/dJ < 0, implying not only an abrupt jump 
but also a hysteresis in E{J), as shown experimentally in 
Ref.1 As B is increased the jump disappears and E{ J) 
is turned to a smooth nonhysteretic curve. 

Previous analyses of the mechanisms that cause the 
FFI relied mostly on identification of the jump at Ei{Ji) 
and discussion of the {B,T) dependences of Ei, Ji and 
other related parameters {ui, power density JiEi, etc.) 
The quantitative description of E[J) over a wide E- 
range, i.e., extending both below and above Ei, has re- 
mained an open question. Furthermore, Ei is a relevant 
parameter even if there is no jump but it has not been ex- 



tracted from the nonhysteretic E{J). The lack of such an 
investigation in conventional superconductors could pos- 
sibly be explained by the usually strong pinning, which 
complicates treatments of pure flux flow effects even in 
simple vortex systems. In high- Tc superconductors the 
pinning is weak, but the exact form of vortices is in this 
case less well known because of their complex anisotropic 
character and peculiar fluctuation phenomena in the de- 
pinned state. ^ To avoid the mentioned obstacles as much 
as possible we have chosen a material already proven to 
be appropriate for studying the fundamental mechanisms 
of vortex dynamics, namely amorphous Nbo.vGeg.a thin 
film of thickness comparable to the coherence length ^ 
These samples have very weak or negligible pinning over a 
considerable part of the {B, T) plane and represent a sim- 
ple classical "dirty" superconductor with a well defined 
vortex structure. In addition, we have reduced the mea- 
surement current and thus the power dissipation in the 
sample by patterning strips of a few micrometre width. 

Close to Tc we have found a convincing quantitative 
agreement with the LO theory up to an unexpectedly 
high b = B/Bc2 ~ 0.7, in both the close-to-equilibrium 
fiux fiow resistivity pf and the E{J) extended over a 
wide range of J. In particular, we show that Ei can 
be determined unambiguously even if there is no clearly 
defined jump. At low temperatures the LO description 
breaks down, which suggests a different physics of the 
FFI. These data can be explained consistently by elec- 
tron heating to a temperature T* above the bath tem- 
perature Tq. The electron heating causes a decrease of 
Bc2{T*) and eventually a transition to the normal state 
at a certain electric field Ec- 

The methods of sample fabrication and determina- 
tion of superconducting parameters are described in 
Ref. [13 The measured microbridge,*^ deposited onto 
an oxidised Si substrate, is 210 /im long, 5 /im wide, 20 
nm thick, and has the following parameters of interest: 
Tc = 2.75 K (the transition width of 0.05 K), the esti- 
mated T — normal state resistivity pn(0) = 3.3 ± 0.2 
^17m, -{dBc2/dT)T=T, « 2.6 TK"*, ^(0) = 6.8 nm, 
and the other Ginzburg-Landau parameters are k = 
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FIG. 1: E{J) at To = 2.5 K, for 0.1 T < B < 0.5 T 
{Bc2 ~ 0.65 ± 0.03 T) increasing as indicated by the arrow, 
measured by sweeping the applied current up. The dashed 
line represents E„ — p„J. Inset: Measured magnetoresistiv- 
ity (solid line) and the LO p/ (open symbols) plotted using 
gi and 172 as explained in the text. 



103 and A(0) = 1.15 fim. All the sample parame- 
ters are within the range of expected values for amor- 
phous Nbo.vGeo.a thin films. The measurements were 
performed in a ^He cryostat with rf filtered leads. The 
dc E{J) was measured by sweeping an applied current 
at a rate 10 nAs~^ (0.1 MAm~^s~^), whereas the mag- 
netoresistivity [p{B, T)] measurements were carried out 
using small currents (1 MA/m^) at which the E{J) is 
linear, originating from either thermally activated vortex 
hopping or free flux flow^ 

In Fig. 1 we show the E{J) at Tq = 2.5 K (t = Tq/Tc = 
0.91) for 0.1 T < S < 0.5 T ( 0.15 <b< 0.77). AU the 
curves were measured by increasing the applied current. 
A change from an E{J) with the FFI jump (low B) to 
a smooth E{J) (high B) is clearly visible, as well as a 
gradual approaching the normal state electric field En — 
PnJ (dashed line) at large J. We show below that the 
LO FFI theory explains quantitatively all these curves. 
Close to Tc the LO expression for J{E) is given by 
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where cr„ — l/pm ^ is a constant of order unity, Y{E) = 
1/(1-1- E'^/Ef) describes the vortex core shrinking, and 
g{b) is a function approximated by the following inter- 
polation formulae: gi{b) = 4.04 - b^/'^{3M + 2.386) for 
b < 0.315, and 32(6) = 0.43(1 - 6)^/2 + 0.69(1 - 6)^/^ for 
b > 0.315. In the limit E < E„ Y{E) w 1 and Eq. [H 
gives the flux flow resistivity pf = E/ J. In the expression 
for Pf, A — I follows from the condition pf{Bc2) = Pn, 
whereas in nonequilibrium the constant value of A w 1 
reflects suppression of the superconducting order param- 
eter outside the cores by a strong electric field.- 

A comparison of the measured data (full lines) and 
Eq. n (dashed lines) is shown in Fig. 2 for two char- 
acteristic shapes of the E{J) at Tq = 2.5 K, i.e., for 
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FIG. 2: E{J) at To = 2.5 K (full lines) measured at an in- 
creasing applied current, for B — 0.1 T (a) and B = 0.4 T 
(b). The dashed lines are plots of Eq.Qwith the appropriate 
choices of the parameters, as discussed in the text. Inset to 
(a): Extracted Ei{B) (circles), illustrating the validity of the 
LO theory with a critical velocity Ui = Ei/B independent of 
magnetic field. 



those with (B = 0.1 T) and without {B = 0.4 T) the 
jump. Eq. n agrees with the data excellently by tak- 
ing (T„ = 3.1 X 10^ S/m from p{B ~ 2Bc2), A ranging 
from 0.92 to 0.97 with no systematic i?-dependence, and 
using Bc2 = 0.65±0.03 T to calculate g{b)/b and the cor- 
responding error bars (important at low b where g{b)/b 
is a steep function). Thus, the only free parameter left is 
Ei, shown in the inset to Fig. 2(a) and discussed later. 

In the inset to Fig. 1 we show p{B) (solid line) at 
the same temperature, compared with the theoretical pf 
(open symbols) . With the same values of parameters i?c2 , 
Gn and A as above, the agreement of the data and the 
LO theory is satisfactory below ~ 0.5 T all the way down 
to _B — > 0. This implies a negligible critical current den- 
sity Jc and a good description of the close-to-equilibrium 
transport properties in terms of the LO theory for all the 
E{J) shown in Fig. 1. The LO theory however fails to 
explain the data closer to i?c2, in contrast to our previous 
flnding^° for another sample at i = 0.82. To resolve this 
apparent inconsistency wc measured p{B) of the present 
sample at Tg = 2 K (t = 0.7) and found that the agree- 
ment of the data close to Bc2 and the LO theory is re- 
stored at this temperature. The failure of the LO theory 
to describe p[B — s- Bc2) in the vicinity of Tc may be re- 
lated to a widening of the equilibrium critical-fluctuation 
region at B sufficiently close to Bc2- 

From the slope of linear Ei{B) we calculate the critical 
vortex velocity Ui = 305 m/s, the quasiparticle energy 
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relaxation timsi t, = D[14C(3)(1 - t)fl'^l'Ku\ = 0.18 
ns near Tf,; where C is the Riemann zeta function, D — 
9.kBTc(^{Q)/TTTi = 4.3 x 10"^ mVs the diffusion constant, 
and the corresponding inelastic relaxation length — 
^/Dtc = 87 nm. The linearity of Ei{B) provides strong 
evidence for the FFI being caused by the LO mechanism 
of vortex core shrinking. Note that the LO model holds 
up to an unexpectedly high h, which is twice larger than 
the upper limit estimated by BS. Only for B — 0.5 T the 
relatively large error bar of the corresponding may 
imply that the BS heating is starting to take place, but 
even there the agreement with Eq.^is very good over the 
whole iJ-range. Previously we have shown that the weak 
heating effects in this regime have contributed mostly to 
the vortex motion noise.«iSi In conclusion to this part, our 
results for Tq close to Tc are over a large _B-interval in 
remarkable quantitative agreement with the LO theory. 

We now turn to the low temperature regime. Recently 
Kunchur analysed the FFI in YBa2Cu307_5 at low tem- 
peratures and small-to-moderate h in terms of electron 
heating to a temperature T*{E) > Tq- The heating oc- 
curs due to insufficient heat transfer to phonons at low 
temperatures, which reduces the superconducting order 
parameter and results in a decrease of Bc2{T*)i' Having 
found the above excellent agreement with the LO pic- 
ture close to Tc, this urged us to carry out Iow-Tq mea- 
surements over a similar range of b and thus investigate 
the differences and/or similarities between the FFI and 
overall nonlinearities of the E{J) at low an high tem- 
peratures. The resulting set of E{J) measurements at 
To = 1.1 K (t = 0.4) is shown in Fig. 3. To check whether 
the LO approach could be valid thus below Tc we tried 
to modify Eq. 1 by using the given form of Y{E) and re- 
placing the ^-dependent part with the ones appropriate 
at low temperatures (see Eas. l2l3l below'). Despite the ap- 
parent similarity of the curves when compared to those of 
Fig. 2, we did not obtain any satisfactory agreement even 
if all the numerical parameters were left floating. This 
motivated us to analyse these results in terms of electron 
heating as the cause of a second type of the FFI. 

At low t and h the expression for J{E) is given bj*i 

Eq. 121 predicts that if the heating is weak the slope of 
E{J) sufficiently above Jc is 0.9/9„6. The dotted lines in 
Fig. 3 show plots of Eq. H for 0.6 T and 0.9 T, where 
the linear region can be found with reasonable reliabil- 
ity. We used Bc2 = 3.0±0.1 T as determined from 
Pf{B) and Jc being the only free parameter. A replace- 
ment To T*(E) in Eq. [21can explain the rise of E( J) 
above the dotted lines by a progressive decrease of Jc 
and increase of b. Moreover, the ratio E/h{E) may de- 
pend nonmonotonically on E. This results in a negative 
slope of J{E) and thus causes another type of flux flow 
instability. However, Eq.|21does not describe the E{J) at 
large E, since the upper limit of the dissipation (given by 
h = 1) cannot explain the data close to This topic 




FIG. 3: (a) E{J) at To = LI K (solid lines), for 0.6 T < B < 
2.2 T {Bc2 = 3.0 ±0.1 T) increasing as indicated by the 
arrow. The measurements correspond to an ascending applied 
current. The dashed line shows En = pnJ- The dotted lines 
are plots of Eq.Hfor 0.6 T and 0.9 T. (b) E vs. b* calculated 
from the measured E[J) and Eq. |H]using a = 3. The vertical 
scale is the same as in (a) and the arrow points again in the 
direction of increasing B. 

was not addressed by Kunchur, who concentrated only 
on the low-5 and \ow-E regime. 

In order to analyse the E{J) En data we recall 
another LO result, namely that as long as the electron 
mean free path is much smaller than ^, close to Bc2 the 
J{E) is at an arbitrary temperature determined bjsi 

J = (7„[l + a(T) (1-5)]S , (3) 

where a is a temperature-dependent constant varying 
between 2 and 4, and Jc at such high dissipation can 
be disregarded. The above expression was used success- 
fully in a detailed analysis of the high-6 magnetoresis- 
tance of amorphous Nbo.yGeo^s filmsii^ If the assump- 
tion of electron-heating-induced nonlinearities is correct, 
Eq. O should describe the upper part of E{J) through 
i?-dependence of b and a up to the transition to the nor- 
mal state at a critical electric field Ec{B) correspond- 
ing to T* = Tc{B) (equivalently, to B = Bc2{T*)). In 
Eq. Othe heating affects a and increases 6 to a nonequi- 
lihrium b*\T*{E)]. The temperature dependence (and 
hence the S-dependence) of a is weakfi and if we approx- 
imate a by a constant we can invert Eq. |31 to calculate 
b*{E) = B/Bc2{T*) from our E[J) data. In Fig. 3(b) 
we show a plot of E vs. b* (calculated using a = 3) for 
h* > 0.9, where we expect the validity of Eq. Oand the 
approximation of a constant a. 

Since the heating characteristics T*(E) is mainly gov- 
erned by electron-phonon scattering it should be inde- 
pendent of B, the effect of which is expressed through 
the Ec{B) dependence. Thus, we expect a scaling be- 
haviour of J{E) for different magnetic fields provided 
Ec{B) is taken into account. This is demonstrated in 
Fig. 4, where Ec{B) is chosen such that b* (calculated 
with a = 3) scales with 1 — E / Ec{B), i.e., the data shown 
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FIG. 4: Scaling plot of the nonequilibrium reduced magnetic 
field b* vs. 1 — E/Ec calculated from Eq.|3]using a — 3, where 
Ec{B) is plotted against equilibrium (1 — b) in the inset. The 
error bars of Ec indicate how much Ec varies if the scaling is 
performed using a between 2 and 4. 

in Fig. 3(b) can be collapsed onto the same curve. The 
corresponding electron temperature isl.7K<r*<2.5 
K, as estimated from the equilibrium Bc2{T) characteris- 
tics. In the inset to Fig. 4 we plot Ec against equilibrium 
1 — 6. The error bars indicate the range of variations 
of Ec when a changes between 2 to 4 (and results in a 
satisfactory scaling). 

The above procedure corresponds to a determination 
of Ec{B). The solid line in the inset to Fig. 4 represents 
Ec = £'co(l — b), where Eco = 1500 V/m, and describes 
the inferred values of Ec fairly well. This Ec{h) depen- 
dence can be related to the thermodynamic properties 
of the mixed state. For a "dirty" high-K superconduc- 
tor the Gibbs free energy density stored in the superfluid 
at sufficiently large h is given by = ?7c(l — &)^, where 
Uc = Bc2/4:fj,oK^ is the zeio-B superconducting condensa- 
tion energypi^In the present case, C/c( 1.1 K) « 170 J/m^. 
Assuming that electron heating weakens the superfluid 
we can equate Gg = anE^Te, which explains Ec <x. 1 — b 



provided does not depend on B (expected for electron- 
phonon scattering). Furthermore, if is considered as 
fairly a constant in the given T*-range we can calculate 
Te(T Tc) ~ Uc/cFnE'^Q — 0.24 US, which is in reasonable 
agreement with the independent estimate = 0.18 ns 
from the E{J) at Tq = 2.5 K. 

Our and Kunchur's complementary approaches agree 
in the final conclusion that electron heating and the 
corresponding reduction of Bc2 is the main source of 
nonlinearities of E{J) in strongly nonequilibrium vor- 
tex transport at low temperatures. Another implica- 
tion of our analysis is that critical fiuctuations around 
Bc2{Ta — > Tc), which affect p/ in complete thermal equi- 
librium (see the inset to Fig. 1), seem to be of little rele- 
vance when To ^ Tc even if i? — s- Bc2 by electron heating. 
This issue however requires a more detailed investigation. 

In conclusion, we have measured and analysed the 
E{ J) curves of amorphous Nbo.TGco.a microbridges over 
a wide range of magnetic field and in two character- 
istic regimes, i.e., close to and well below Tc. In the 
former case we have found an excellent agreement with 
the Larkin-Ovchinnikov theory of nonlinear flux flow and 
the related instability up to a surprisingly high value of 
B/Bc2 ^ 0.7, much larger than predicted theoretically 
by Bezuglyj and Shklovskiij. At low temperatures the 
nonlinear! ty oi E{J) and the flux flow instability can be 
reasonably well described by electron heating and the re- 
lated decrease of Bc2- Our scaling analysis of the E{J) 
curves at high currents supports this conclusion quantita- 
tively, leading to an agreement with the thermodynamic 
properties of the mixed state. 
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